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We report on a theoretical study of magnetic transitions induced by tunnelling electrons in 
individual adsorbed M-Phthalocyanine (M-Pc) molecules where M is a metal atom: Fe-Pc on a 
Cu(110)(2x l)-0 surface and Co-Pc layers on Pb(lll) islands. The magnetic transitions correspond 
to the change of orientation of the spin angular momentum of the metal ion with respect to the 
surroundings and possibly an applied magnetic field. The adsorbed Fe-Pc system is studied with 
a Density Functional Theory (DFT) transport approach showing that i) the magnetic structure 
of the Fe atom in the adsorbed Fe-Pc is quite different from that of the free Fe atom or of other 
adsorbed Fe systems and ii) that injection of electrons (holes) into the Fe atom in the adsorbed 
Fe-Pc molecule dominantly involves the Fe 3d z 2 orbital. These results fully specify the magnetic 
structure of the system and the process responsible for magnetic transitions. The dynamics of the 
magnetic transitions induced by tunnelling electrons is treated in a strong-coupling approach. The 
Fe-Pc treatment is extended to the Co-Pc case. The present calculations accurately reproduce the 
strength of the magnetic transitions as observed by magnetic IETS (Inelastic Electron Tunnelling 
Spectroscopy) experiments; in particular, the dominance of the inelastic current in the conduction 
of the adsorbed M-Pc molecule is accounted for. 



PACS numbers: 68.37.Ef, 72.10.-d, 73.23.-b, 72.25.-b 

I. INTRODUCTION 

The development of Scanning Tunnelling Microscopy 
(STM) recently enabled to study the case of localised 
magnetic excitations at surfaces 1 - - — . Low-temperature 
STM experimental studies revealed that, in certain cases, 
a local spin could be associated with individual adsor- 
bates at surfaces. Under the variation of the STM bias, 
the tip-adsorbate junction exhibits conductance steps at 
well-defined energies in the few meV energy range. These 
steps are attributed to the excitation of the local spin of 
the adsorbate and the step energy position yields the cor- 
responding excitation energy, leading to magnetic IETS 
(Inelastic Electron Tunnelling Spectroscopy). Reveal- 
ing the existence of nano-magnets at the atomic (molec- 
ular) level can have fascinating consequences for the 
miniaturisation of electronic devices. In addition, the 
possibility to determine the energy spectrum of the lo- 
cal spin on the adsorbate as a function of an applied 
magnetic B field provides an efficient way to quanti- 
tatively characterize the local spin and its interaction 
with the underlying substrate. This makes magnetic 
IETS an invaluable tool for magnetic studies of nano- 
objects at surfaces. In this respect, one can mention 
several spectacular results obtained in this way: evi- 
dence of spin coupling between neighbouring atomic ad- 
sorbatcs (anti-ferromagnetic coupling along adsorbed Mn 



chains)^, large magnetic anisotropy of atomic adsorbates 
(Mn and Fe adsorbates on CuN) 3 -, change of magnetic 
structure of a molecular adsorbate on various substrates 
(Fe-Phthalocyanine on Cu and CuO)^, evidence of super- 
exchange interactions^, interaction of molecular adsor- 
bates with magnetic substrates^, charging of adsorbed 
magnetic nano-objects 7 -. Magnetic transitions could be 
easily evidenced in these systems due to the presence of 
a coating on the metal substrate that efficiently decou- 
ples the adsorbate from the continua of metallic states. 
When magnetic atoms are directly adsorbed on a metal, 
magnetic transitions could be observed but much broad- 
ened by the interaction with the substrate^. 

Besides detailed insight into the structure of a mag- 
netic adsorbate on a surface, IETS experiments bring also 
information on the dynamics of the excitation of a local 
spin by tunnelling electrons. Indeed, when considering a 
conductance spectrum as a function of the junction bias, 
the position of the conductance steps yields the excited 
state energies and the height of the conductance steps 
relative to each other and relative to the conductance at 
zero bias yields the relative magnitude of the various pos- 
sible magnetic excitations. This feature is quite impor- 
tant for possible future applications, since it determines 
how easily a local adsorbate spin can be flipped at will 
by tunnelling electrons or can be quenched by collisions 
with substrate electrons^. It appears that magnetic tran- 
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FIG. 1: Scheme of an M-Phthalocyanine molecule. In the 
present study M is taken as a Fe or a Co atom, here repre- 
sented in red. The free molecule has a D±h symmetry which 
is generally reduced upon adsorption on a substrate. The N- 
atoms are depicted in blue, C atoms in grey and H atoms in 
cyan. 

sitions are highly probable^—, much more probable than 
other inelastic processes like e.g. vibrational excitation 
by tunnelling electrons^— . As an extreme example, 
in the case of Co-Phthalocyanine (Co-Pc)£, the inelas- 
tic contribution to the current is found to be three times 
larger that the elastic current, whereas for vibrational 
excitation, the inelastic tunnelling contribution reaches 
at most a few per cen^r— . 

A few tools have been used to describe the spec- 
troscopy and dynamics of a local spin on an adsorbate. 
First, the energy of a local spin interacting with its envi- 
ronment and possibly with an applied magnetic field, B, 
has been modelled efficiently with the following magnetic 
Hamiltoniani^: 



introducing in the description of electron tunnelling a 
coupling term proportional to the local spin. In first 
order perturbation theory, the inelastic current then ap- 
pears proportional to the squared matrix elements of the 
coupling term between initial and final states, similarly 
to what was noticed in Ref. [||. Recently, a completely 
different approach was introduced to treat magnetic exci- 
tations by tunnelling electrons^. It is based on the large 
difference of time scale between electron tunnelling and 
magnetic anisotropy of adsorbates: electron tunnelling 
is fast and the magnetic anisotropy can be considered 
as non-active during tunnelling. One can then describe 
tunnelling without the magnetic anisotropy taken into 
account and simply switch the latter at the beginning 
and at the end of tunnelling. In addition, DFT calcula- 
tions on the two studied systems (Fe and Mn adsorbates 
on CuN) showed that only one coupling scheme between 
adsorbate and tunnelling electron spins is significantly 
contributing to tunnelling. The magnetic excitation ap- 
pears as the result of a spin decoupling/rccoupling pro- 
cess. This approach is non-perturbative and can then 
handle the large excitation probabilities encountered in 
these systems as well as make predictions for the elas- 
tic/inelastic relative contributions to tunnelling. In ad- 
dition to bringing a qualitative view into the magnetic 
excitation process, it was shown to accurately account 
for the strength of the magnetic excitations, in particu- 
lar relative to the elastic channel, in the case of Mn and 
Fe atomic adsorbates on CuNil. In the present work, we 
show how our method can be used in the case of Fc-Pc 
and Co-Pc molecules (Fig. [T]) adsorbed on partly insulat- 
ing substrates. As one of the main results, this approach 
is shown to precisely account for the overwhelming dom- 
inance of the inelastic conductance over the elastic con- 
ductance observed experimentally in these systems^. 



H = gfi B B-S + DS 2 + E(S 2 ~S 2 ), (I) 

where S is the local spin of the adsorbate, g the Landc 
factor and \ib the Bohr magneton. B is an applied mag- 
netic field. D and E are two energy constants describing 
the interaction of S with the substrate, i.e. the mag- 
netic anisotropy of the system. The Cartesian axis xyz 
are chosen according to the magnetic symmetries of the 
system. Hamiltonian (fTJ) describes how the adsorbate 
spin, S, is oriented in space due to its interaction with 
the substrate. The observed magnetic excitation energies 
of adsorbates were very efficiently modelled by the above 
Hamiltonian^—. As for the strength of the magnetic exci- 
tations, it was first analyzed in a phcnomcnological way 
by Hirjebehidin et afi. They showed that the relative 
heights of the inelastic conductance steps were very close 
to the relative magnitude of the squared matrix elements 
of the operator between the initial and final states of the 
transition. Nothing could be said about the relative mag- 
nitude of the elastic and inelastic conductance. This find- 
ing was later supported by several theoretical studies^— 



II. METHOD 

We consider a local spin, S, localized on an adsor- 
bate on a surface. It is coupled to its environment 
by the Hamiltonian (JfJ. Diagonalizing the Hamilto- 
nian (fTJ) yields the various spin states of the system, 
their energies, E n , and the associated wave-functions, 
(f) n . Hamiltonian (fTJ) is more easily written in the ba- 
sis of \S, M) states, the eigenstates of S 2 and S z , so that 
the anisotropy states can be written as: 

\<t>n) =J2 C n,M\S,M) (2) 
M 

In the present study, we used the D, E and g parameters 
in Hamiltonian (JXJ) as determined in the analysis of the 
experimental energy spectrum and from these, without 
any further parameter adjustment, we derive the strength 
of the magnetic transitions induced by tunnelling elec- 
trons. The various energy terms in Hamiltonian (fTJ) arc 
in the few meV range, so that one does not expect the 
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FIG. 2: Projected density of states (PDOS) on the Fe d- 
atomic orbitals. For all the curves shown here, the positive 
(black) curves corresponds to the majority spin, and the neg- 
ative (red) to the minority spin. The d orbitals are classified 
according to the cartesian axes that contain the N-Cu-N axis 
of the molecule (x,y,z), or with respect to the surface direc- 
tions: x' for the [1 — 10] and y' for the [001] directions. The 
z axis is the same for both reference frames. 



corresponding interaction to play a role during the elec- 
tron tunnelling process and one can treat this problem 
in the sudden approximation: treat the tunnelling with- 
out the magnetic anisotropy and then introduce it as a 
frame transformation at the beginning and at the end of 
tunnelling. 

Electron tunnelling from the STM tip through the ad- 
sorbatc and into the substrate (in the absence of magnetic 
anisotropy) can be represented by a scattering T matrix, 
noted Txip^Sub (and an equivalent one for the reverse 
tunnelling). It corresponds to the scattering of the tun- 
nelling electron by the adsorbate and it depends on the 
electron energy. In the absence of magnetic anisotropy, 
it depends a priori on the spin coupling between the tun- 
nelling electron and the adsorbate, via the exchange in- 
teraction. Thus, in the absence of significant spin-orbit 
interactions, it can be written in a diagonal form if we 
consider St, the total spin of the system (electron + ad- 
sorbate). Defining \St,Mt) as the cigenfunctions of S T 
and St,z (if S is the adsorbate spin, then St = S ± |), 
we can write formally the scattering TTip^sub matrix (in 
the absence of magnetic anisotropy) as: 

T T i P ->sub = \^T^ M T)T T T p _^ Sub (S T ,M T \. (3) 

S T ,M T 



is a number, function of the electron energy. 
However, since, below, we consider only the limited en- 
ergy range spanned by the magnetic excitations (up to 
10 meV), the transmission probabilities [T^lp^su^ 2 can 
be considered as constant in the present study. In the 
sudden approximation, the tunnelling amplitude (in the 
presence of magnetic anisotropy) is written as the matrix 
element of the TTip^Sub amplitude between the initial 
and final states of the tunnelling process. These states 
are written as |i,m;0 n ) where the first part concerns 
the tunnelling electron (the electron spin is 1/2 and m 
is the projection of the tunnelling electron spin on the 
quantization axis) and the second part concerns the local 
spin of the adsorbate. One then obtains the amplitude, 
AM P mi n-^m' ,n' , for a tunnelling electron induced tran- 
sition from (f>„ to <f) n i , while the tunnelling electron spin 
projection changes from m to ml as: 



AMP rn 

X TjM t {\ 



St 



m';<p' n \Srj 



— J2s T ^-Tip- 

,M t )(St,Mt\ 



Sub 

hm;<f> n ) 



(4) 



One can see that there is a superposition (interference) 
of tunnelling amplitudes through the different St intro- 
duced by the magnetic anisotropy. Equation yields 
the transition amplitudes in the general case, if the spin 
of the tunnelling electron is registered in both the initial 
and final states. Implicitly, it has been assumed above 
that the tunnelling electron quantization axis is the z-axis 
of the adsorbate magnetic anisotropy; situations with dif- 
ferent quantization axis for the adsorbate and the tun- 
nelling electron can be easily handled with an expression 
similar to Eq. 

We can now define the probability, P, for transitions 
from <f) n to (j> n > induced by unpolarized tunnelling elec- 
trons by summing incoherently over the distinguishable 
channels: 



P, 



m,m' 

771,771' St 

x £(1 ml- 4>' n \S T , M T )(S T , M T \\,m; <j> n )\ 2 . 



Mr 



We can be a little more explicit by writing the total 
spin states, \St, Mt), as expansions over uncoupled spin 
states: 



| St, Mr) = J2 CG s T ,M T , m \S, M 



Mt 



■ m )\^, m ) (6) 



where |S, M) states correspond to the adsorbate spin 
states and | ^ , m) to the tunnelling electron spin. Again, 
the adsorbate and electron spins are quantized on the 
same axis. The CG are Clebsch-Gordan coefficients. 
Combining ^ and © we can express the total spin 
states as an expansion over uncoupled products of adsor- 
bate magnetic anisotropy states and tunnelling electron 
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spin: 



\3) 



\S T ,M T ) 



\4> n )\\/2,m) (7) 



Equation (0 links the \j) = \St, Mt), states appropri- 
ate to describe tunnelling without magnetic anisotropy to 
the channel states of the complete tunnelling process. We 
can then rewrite the transition probability (j6j as: 



Pn—yn' — 



ErpS T 
± Tip^Sub 

St 



(8) 



In the case where the T T J p ^ Sub tunnelling amplitude is 
dominated by a single symmetry, St (like it was found 
in the case of Mn and Fe adsorbates on CuN — , and like 
it is shown below to be the case in the systems studied 
here), the probability (JSJ further simplifies into: 



P, 



I rpST 

2 l^Tip^Subl 



2 E 



M T 



This result, used in Ref. [17j . is very simple, the elec- 
tronic part of the tunnelling (the T T J 3ub amplitude) is 
factored out and the probabilities for the different chan- 
nels are simply proportional to spin-coupling coefficients 
corresponding cither to the magnetic anisotropy or to the 
coupling between electron and adsorbate spins (the co- 
efficients are products of the diagonalization expansion 
coefficients in Eq. ([2]) and Clcbsch-Gordan coefficients). 

From the transition probabilities we can write the con- 
ductance dl/dV as a function of the STM bias, V, as: 



dl_ 

dV 



Co 



£„ Q(V - EX n ) £ mim , | A jXn Al njn , | 2 



V A ■ , A* 2 



(10) 

Expression (|T0|) corresponds to the conductance for the 
system being initially in the ground state n = 1. The 
sum over n extends over all the \4> n ) states, including the 
ground state, so that the above conductance takes all con- 
tributions, elastic and inelastic, into account. EX n is the 
excitation energy of the magnetic level n, corresponding 
to the eigenvalue difference of the final, \<j> n ), and ini- 
tial \<pi) states. The Heavyside function, G, takes care of 
the opening of the inelastic channels at zero temperature. 
Co is the total conductance corresponding to the trans- 
mission amplitude T^f _^ Sub . It is then a magnetism- 
independent conductance. Since we only consider a lim- 
ited V range, defined by the magnetic excitation energies, 
Cq can be considered as constant in the relevant V-rangc. 
Co is equal to the conductance of the system for biases 
larger than all the inelastic thresholds. Expression (|T0|) 
corresponds to the case where only one St value actu- 
ally contributes to tunnelling so that the sum over j is 
restricted to the corresponding Mt values. If the two 
St symmetries contribute to tunnelling, a more general 
expression derived from Eq. ((SJ has to be used. 

The above form © for the transition probabilities ac- 
counts for the strength of the magnetic transitions. In 



fact, the transition probability is not proportional to a 
coupling term as in perturbation approaches, it appears 
as a simple sharing of the global transmission of the tip- 
substrate junction, \T^l p ^ Sub \ 2 , among the various spin 
states, \<fi n ). Qualitatively, in the magnetic excitation 
process, one can say that the incident electron arrives in 
a given spin state, it couples with the adsorbate spin to 
form a state of the total spin, St', tunnelling through the 
junction occurs independently in the different total spins; 
at the end of the electron-adsorbate collision, the total 
spin splits back into its adsorbate and electron compo- 
nents, populating all the possible adsorbate spin states. 
Then, the expression for the sharing process (jH]) simply 
expresses angular-momentum conservation. 



This kind of excitation process has been invoked in 
various situations where an exchange of angular momen- 
tum is involved: resonant rotational excitation in elec- 
tron collisions on free and adsorbed molecules 18 i 19 , spin- 
forbidden transitions in electron-molecule collisions 2 ^ or 
in atom-surface scattering 2 ^. In all cases, it leads to ef- 
ficient excitation processes. One can stress that this de- 
scription of magnetic transitions is at variance with that 
of the vibrational excitation induced by collisional elec- 
trons— ^ or by tunnelling electrons 2 ^—. However, even 
in the vibrational excitation case, a process similar to 
the one discussed here for magnetic transitions does ex- 
ist. It is associated with the composition of the incident 
electron linear momentum with that of the target atom. 
Conservation of momentum leads to recoil of the target 
induced by electron scattering; however, due to the large 
ratio between electron and nuclei masses, the correspond- 
ing excitation is very weak and leads to very small exci- 
tation probabilities. Vibrational excitation has then to 
involve other processes such as e.g. resonant scattering, 
where the increase of the collision time allows weak in- 
teractions to be efficient in various situations of electron- 
molecule collision o 22 ' 23 ' 27 " — . The above discussion can 
be summarized by stressing that the angular momenta 
(orbital or spin) of electrons and atoms or molecules can 
be of the same order of magnitude due to quantization 
(a few units in many cases) so that exchange between 
them is easy and has visible effects; in contrast, because 
of the large electron-nuclei mass ratio, the electron linear 
momentum is usually much smaller than that of heavy 
particles limiting the efficiency of momentum exchange 
processes. As a consequence, angular and spin degrees 
of freedom appear similar and the present treatment of 
spin transitions is very similar to the treatment of ro- 
tational excitation used in Ref. [f9[ to account for the 
experimentally observed efficiency of tunnelling electron 
in inducing molecular adsorbate rotations^ 2 -. 
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FIG. 3: Amplitude of the Kohn-Sham orbital of the Fe- 
Pc/Cu(110)(2xl)-O + STM tip system at 0.2 eV above the 
Fermi level (positive in light grey and negative in light pink). 
The STM tip in the upper part is made of a Cu(110) surface 
with an extra protruding atom (in gold) . The adsorbed Fe-Pc 
is lying flat on the substrate (N atoms in blue, C atoms in 
grey and H atoms in cyan). The substrate is in the lower part 
of the figure (Cu atoms in gold and O in red). The orbital is 
concentrated around the Fe atom and exhibits a strong d z 2 
character perturbed by the interaction with the substrate. 



III. MAGNETIC EXCITATIONS IN 
SUPPORTED Fe-PHTHALOCYANINE 
MOLECULES 



A. Electronic structure of the Fe-phthalocyanine 
molecule 

A detailed description of the electronic structure of the 
Fe-Pc can be found in various references (see e.g. [33l - [35j 
and earlier references there in). The s outer electrons 
of the Fe atom are transferred to the surrounding Pc 
ring leaving a central Fe 2+ ion with a 3d 6 electronic con- 
figuration. In the free Fe-Pc, see Fig. [1] the Fe atom 
is surrounded by the Pc ring of symmetry, so that 
the Fe d manifold splits into &2g (d xy ), e g (d zx ,d yz ), a\ g 
(d z 2) and b\ g (d x 2_ y 2) orbitals (the z-axis is normal to 
the Pc plane and the x and y-axis are parallel to the 
Fe-N interatomic axis ). In the free Fe-Pc, the bi g {d xy ) 
orbital is fully occupied, the b\ g (d x 2_ y 2) corresponds to 
the highest eigenenergy because of the large overlap with 
the N-atom orbitals, and four electrons occupy the e g 
(dzxjdyz) and a\ g (d z 2) orbitals. Various configurations 
have been proposed, with small energy differences be- 
tween them (see a discussion in Ref. |35(). The structure 
of the other M-Pc (M= metal) is similar—, a change of M 
along the Fe, Co, Ni, Cu and Zn sequence corresponding 
to the further filling of the split d manifold. 

When the Fe-Pc is adsorbed on a Cu(110)(2x l)-0 sur- 
face, the symmetry is broken. If we assume the pres- 
ence of the surface to only induce a perturbation on the 
structure of the free Fe-Pc, then the g subscript of the 
orbitals disappears and the e g orbital splits (the d zx and 
d yz are not degenerate anymore) . We will use the orbital 
notation of the free Fe-Pc when discussing the electronic 
structure of adsorbed Fe-Pc, although the molecules are 
distorted by the interaction with the surface. 



Recently Tsukahara et a& performed a detailed mag- 
netic IETS study of single Fc-Phthalocyanine (Fe-Pc) 
molecules adsorbed on a Cu(110)(2xl)-O surface. The 
molecule lay flat on the surface and two adsorption ge- 
ometries were found, labelled a and /3 differing by the rel- 
ative orientation of the molecule on the substrate. Clear 
magnetic transitions were observed by scanning the STM 
tip bias and were attributed to a local spin S = 1 inter- 
acting with the environment and with an applied mag- 
netic field, B. The magnetic transitions were only ob- 
served when the tip was placed above the Fe atom and 
were attributed to a local spin of the Fe atom. The en- 
ergies of the magnetic levels of the system as a func- 
tion of the applied field were very precisely accounted 
for in Ref. ||| using Hamiltonian ([T]). The parameters of 
the two, a and /? adsorption geometries are different: 
D = -3.8 mcV, E = 1.0 meV and g = 2.3 for a-Fc-Pc 
and D = -6.9 mcV, E = 2.1 meV and g = 2.4 for (3- 
Fe-Pc, i.e. the same kind of structure but with a very 
different zero-field splitting of the magnetic states. We 
used this modelling of the magnetic structure to com- 
pute the strength of the magnetic excitations as induced 
by tunnelling electrons, making use of the formalism de- 
scribed in section HIl 



B. Density functional study of Fe-Pc adsorbed on a 
Cu(110)(2xl)-O surface 

The ground state electronic structure configuration 
and the global transmission, T^f _^ Sub for the tunnelling 
from an STM tip to the substrate passing through the 
Fe-Pc molecule were obtained by density functional the- 
ory (DFT) simulations, performed with the Siesta and 
Transiesta code a 36 ' 37 . We have used the generalized 
gradient approximation^, for the exchange-correlation 
functional. We only studied the j3 adsorption geome- 
try of the Fe-Pc on Cu(110)(2x l)-0, which has a higher 
symmetry than the a geometry. The electronic structure 
in the a and j3 geometries are assumed to be equivalent. 

Starting with the electronic structure analysis for the 
adsorbed molecule, Siesta uses atomic orbitals for the 
basis set, and the projection of the density of states onto 
the Fe atomic orbitals reveals an open shell structure 
with two unpaired electrons, i.e. a S = 1 configuration 
{dl y dy, z dl 2 dl, z ), as shown in Fig. [5] Note that the 
P adsorption geometry corresponds to the Fe-N axis ro- 
tated by 45° from the symmetry axis of the Cu(110) sur- 
face. As a consequence, the splitting of the e g orbital by 
the substrate field involves x' and j/'-axis rotated by 45° 
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FIG. 4: In panel (a), the transmission function for the system 
shown in Fig. [3] as a function of the tunnelling electron energy 
(referred to the Fermi energy Ef) for a vanishing STM bias. 
Panel (b) shows the PDOS on the Fe d z i atomic orbital. Panel 
(c) shows the projection onto the a\ g molecular orbital of the 
free distorted Fe-Pc that corresponds to the d z 2 Fe orbital. 
Positive (black) curves correspond to the majority spin, and 
the negative (red) ones to the minority spin. 



from the x and y-axis that correspond to the splitting 
of the d manifold by the interaction with the Pc ring. 
The projections in Fig. [2] were then performed on the 
appropriate symmetry orbitals i.e. the d xy , d z 2, d x 2_ y 2 1 
d zx i and d y i z orbitals. More quantitatively, the computed 
spin polarization is 1.85 unpaired electrons. The corre- 
sponding spin is roughly S = 1, in good agreement with 
the experiment^. Due to its large overlap with the sur- 
face electronic structure, the d z 2 orbital is the one that 
hybridizes the most with the substrate. 

Figure [3] shows the real space plot of the Kohn- 
Sham orbital at 0.2 eV above the Fermi level for the 
molecule+surface+STM tip system. This orbital corre- 
sponds to the peak in the PDOS onto the d z 2 Fe orbital 
(see Fig. [2]). The strong d z 2 character of the orbital seen 
in Fig. [3] confirms that this orbital corresponds to the 
dig orbital perturbed by the substrate. The intuitive no- 
tion that, by reaching further out along the z-axis, this 
orbital would contribute more to the junction transmis- 
sion is supported by comparing the transmission func- 
tion (panel (a) of Fig. gj) with the PDOS onto the d z 2 Fc 
atomic orbital (b), and PDOS onto the molecular orbital 
a\ g of the Fe-Pc (c). Note that, for this projection, we 
used the a\ g orbital of the free Fe-Pc molecule, with the 
distorted atomic configuration of the adsorbed molecule. 
From this, it is clear that, at the Fermi level, the trans- 
mission is essentially made of the tail of the a\ g resonance 






FIG. 5: Transmission eigenchannel corresponding to the 
largest transmission amplitude for (a) E — Ef where Ef 
is the system's Fermi energy and (b) E — Ef + 0.2 eV. Light 
grey (light pink) color corresponds to the positive (negative) 
imaginary part of the eigenchannel amplitude coming from 
the STM tip. In gold color, the positive real part. Note that 
the isosurfaces were chosen different in (a) and (b) because 
of the large difference in transmission probability. The trans- 
mission channel exhibits a strong ai g (d, 2 ) character around 
the Fe atom. 



due to the d z 2 Fe orbital, and that only the minority spin 
contributes to the tunneling tranmission. In the trans- 
mission calculations, the tip-molecule distance was set at 
a small distance, 5 A, with a large enough transmission 
probability to visualize the electron transmission eigen- 
channel easily. 

This conclusion is further confirmed by computing the 
transmission eigenchannels. For that, we have used the 
Inelastica code — . The d z 2 orbital is the largest con- 
tributor to S-matrix eigenchannel dominating the trans- 
mission of electrons from the STM tip (here represented 
by an atom on a Cu(110) semi-infinite electrode) placed 
above the molecule. As can be seen in Fig. [5j the 
character of the dominant transmission eigenchannel 
is the same at the Fermi energy, Ep, Fig. Eta) and at 
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FIG. 6: Computed conductance for the (a) a and (b) (3 config- 
uration of an adsorbed Fe-Pc molecule on a Cu(110)(2xl)-O 
surface. The conductance has been normalized to 1 at large 
bias and the curves for the various B fields have been offset 
by 0.5 for representation purposes. The B field is oriented 
vertical to the sample with magnitude 0, 3, 6, 9 and 11 T as 
indicated in the graph. 



E = Ep + 0.2 eV, Fig. [5] (b), the peak of the trans- 
mission resonance. These plots nicely correspond to the 
density plot of the Kohn-Sham eigenstate, Fig. [3] 



C. Orbital picture of the electron tunnelling 
process 

From the above DFT study, we can conclude that 
the Fe-Pc on a Cu(110)(2x l)-0 surface, is associated to 
the triplet electronic configuration d 2 . y d y , z d 1 z2 dl,, z consis- 
tently with experiment — . The interaction between the 
d manifold and the surroundings is seen to split com- 
pletely the d manifold, leading to five non-degenerate or- 
bitals; the orbital angular momentum is then completely 
quenched, in the absence of significant spin-orbit cou- 
plings^ and this justifies the discussion, used in Rcf^ 
and here, of the magnetic anisotropy in terms of the 
spin angular momentum orientation. One can also no- 
tice that the Fe spin is also partially quenched from its 
free atom value (S = 2), again due to the interaction 
with the surroundings that induces a large upward energy 
shift of the d x 2_ y 2 orbital. The present DFT study also 
shows that when an STM tip is placed above the Fe atom, 
tunnelling dominantly involves the d z 2 (ai g ) orbital. So 
when an electron is sent from the tip on the Fe, it involves 
the d 2 xy d y , z d 2 z2 d\., z transient configuration, i.e. the total 
spin of the electron-adsorbate scattering intermediate is 
St = 1/2. Similarly, if a hole is sent from the tip to the 
Fe atom, tunnelling involves the d 2 y d 2 , z d° z2 d\, z transient 
configuration, i.e. the total spin of the electron-adsorbate 
scattering intermediate is again St = 1/2. Thus, of the 
two possible symmetries for the electron-adsorbate scat- 



tering (St = 1/2 or 3/2), St = 1/2 is the prevailing 
symmetry in the tunnelling process. 

In addition, transmission through the d z 2 orbital will 
dominate in a constant-current STM image and generate 
a bright spot at the Fe centre. This is consistent with the 
observation in Ref. Q of the Fe-Pc molecule as a bright 
spot at the centre of a clover leaf. The bright spot corre- 
sponds to the d z 2 Fe orbital and the clover leaf is given by 
the contribution from other orbitals localized on the Pc 
ring (see in [34[ a discussion of the link between bright M 
atoms in M-Pc STM images and their d orbitals). Fur- 
thermore, the magnetic transitions in this system were 
also found^ to be localized in the same region about the 
Fe atom, which we attribute to the d,2 orbital. 



D. Magnetic excitation processes 

The inelastic conductance of the Fe-Pc has been com- 
puted using expression (|10[) for the a and j3 adsorption 
geometries, using the spin parameters determined in the 
DFT study. The conductance is shown in Fig. |H{a) for 
a configuration and Fig. [5] (b) for the (3 configuration, 
as a function of the STM bias for various values of the 
applied magnetic field, £?, along the z-axis. The con- 
ductance has been normalized to 1 at large bias and the 
curves for the various B fields have been offset by 0.5. 
In the calculation, a Gaussian broadening of 0.25 meV 
was introduced to mimic various broadening effects. The 
conductance spectra resemble very much those measured 
by Tsukahara et with well-marked steps at the mag- 
netic excitation thresholds (two inelastic thresholds for 
this S = 1 system because anisotropy splits the triplet 
state in three states^) and very significant contributions 
from the inelastic currents. Indeed, the conductance is 
dominated by inelastic tunnelling at large bias! 

The differences between the a and ft adsorption ge- 
ometries are also well reproduced. One can also stress 
that, in the present approach, the magnetic excitations 
do not modify the maximum value of the global conduc- 
tance Co, which is then independent of changes in the 
magnetic structure and in particular independent of the 
applied B field. This appears clearly in Figs. [6] (a) and 
(b) as well as in Fig. 2a and 2b of Ref. [J], confirming 
our view of the magnetic excitation process as a shar- 
ing of the global conductance over the various possible 
magnetic channels. 

The relative magnitude of the elastic and inelastic 
channels in this system are further illustrated in Fig. [7] 
which presents the magnitude of the inelastic conduc- 
tance steps for Co = 1. The two inelastic step heights 
are a.\ and a 2 (ct\ for the lowest threshold) for the a 
geometry (and similarly for the f3 one). The elastic con- 
ductance is then equal to the global conductance minus 
the inelastic ones, (1 - a± - a?) with this definition. 

The present theoretical results as functions of the ap- 
plied B field are compared with the data of Tsukahara et 
a&. As a first remark, inelastic tunnelling contributes sig- 
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FIG. 7: Relative inelastic step heights in the conductance as 
a function of the magnetic field B, for (a) a and (b) j3 con- 
figurations. Qi and a.2 refer to the first and second excitation 
steps for the a configuration respectively. Analogously, j3\ 
and /?2 refer to the first and second excitation steps for the /3 
configuration. The experimental data points are represented 
with black squares for the first excitation and as red circles for 
the second one and are taken from the supplemental material 
of Ref. n. 



nificantly to the tunnelling current: typically for B = 
T, the inelastic current is equal to twice the elastic cur- 
rent. Second, we can see that the present theoretical 
results reproduce very well the relative magnitude of the 
three contributions to the conductance (elastic and two 
inelastic). In particular, the variation with B is well ac- 
counted for. 



The variation of the strength of the magnetic excita- 
tion with B is in fact reflecting the change of the magnetic 
structure of the adsorbate. For B = 0, the structure of 
the conductance curve is due to the anisotropy imposed 
by the substrate to the Fe-Pc molecule; although the level 
positions arc different in the a and j3 adsorption geome- 
tries, the eigenstates of Hamiltonian ([1]) at B = are 
the same in both geometries (see Tsukahara et a&) and 
consequently, the excitation probabilities are the same, 
as seen in Fig. [7] The effect of a finite B field is to de- 
couple the Fe-Pc spin from the substrate and to tend to a 
Zeeman limit at large B. This decoupling is easier for the 
a-gcometry, because of a weaker transversal anisotropy 
E (given above). In the limit of very large B (not fully 
reached here), the <f) n states reduce to \S, M) states with 
the ground state corresponding to M = —1. The only 
allowed transition is to the (M = 0) excited state, so 
that only one inelastic step remains in the conductance 
spectrum. Its height is given by the modulus square of a 
Clcbsch-Gordan coefficient and in this limit, the inelastic 
current is equal to one third of the total current. 



IV. MAGNETIC TRANSITIONS IN 
SUPPORTED Co-PHTHALOCYANINE 
MOLECULES 

Superposed layers of Co-Pc molecules on Pb(lll) is- 
lands were studied experimentally by magnetic IETS and 
revealed the existence of superexchange interactions^. 
Magnetic excitation of the Co-Pc molecules were ob- 
served but only for several adsorbed layers of Co-Pc 
molecules. No excitation was observed in the single layer 
case, suggesting that either the spin of the Co-Pc molec- 
ular layer lying directly on the Pb surface was quenched 
by the interaction with the substrate or the possible mag- 
netic excitations were very short livedo. The Co-Pc spin 
was found to be S = 1/2, and when several layers are 
stacked one on top of each other, with a stacking angle 
of 60°, the spins of the Co-Pc molecule in the second 
and outer layers couple together in a way well-described 
by a Heisenbcrg Hamiltonian, with an anti-fcrromagnctic 
exchange coupling, J, of the order of 18 meV— . 

Very strong magnetic excitations were reported in the 
multi-layer case^. The present formalism can be used 
to predict the magnitude of the magnetic excitation in 
such systems. If we neglect the magnetic anisotropy of 
adsorbed Co-Pc, the only ingredients in our approach 
are the adsorbate spin S, the intermediate total spin St, 
and the anti-ferromagnctic interaction J for the stacked 
molecules. Both S and J are known from experiment^. 
We did not perform a detailed DFT study as for the 
Fe-Pc case. We assumed that we can extrapolate the 
electronic structure from Fe-Pc to Co-Po^, in both cases 
the Pc molecules being partly decoupled from the under- 
lying metal. As discussed in Ref. [3a, an M-Pc series is 
formed by the various incomplete d-shell metals and the 
Co structure corresponds to adding a d electron to the 
Fe case. With the open shell structure of Fc-Pc outlined 
above {d 2 xy d 2 y , x d^d^.i z ) , this leads to a doublet configu- 
ration of the Co-Pc electronic structure, independently 
of the orbital in which the electron is added. This is 
perfectly consistent with the experimental observation. 
Then, we can assume that tunnelling involves minority 
spin electrons because the molecular electronic structure 
at the Fermi level is given by the last partially-occupied 
d orbital also in this case. Hence, the intermediate spin 
state is St = 0, both for electrons and holes. 

In the case of a single active Co-Pc molecule (two ad- 
sorbed layers) , an applied magnetic field is needed to per- 
form a magnetic IETS experiment, because of the molec- 
ular spin S = 1/2, see Ref. For a finite B field, a Zee- 
man splitting corresponding to a Lande factor 1.88 has 
been observed — . Assuming a simple Zccman structure, 
i.e. no magnetic anisotropy induced by the substrate, the 
above formalism, at finite B, predicts a single inelastic 
peak with an inelastic current equal to the elastic current 
(note that this result is independent of the value of the 
Landc factor). This is consistent with the observations^, 
which reported an increase of the current around 110-120 
% at the inelastic threshold (see Fig. 2c in Ref. (H). 
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In the case of three molecular layers, two stacked 
molecules in the upper two layers interact via an anti- 
ferromagnetic coupling, so that the ground state of the 
two-molecule system is a spin state while the spin 1 
states are excited states (we neglect a possible anisotropy 
induced by the substrate). Magnetic excitations of the 
system by tunnelling electrons for a vanishing B have 
been observed experimentally as a sharp step in the con- 
ductance^. The above formalism predicts at B = 0, in 
the absence of any anisotropy induced by the substrate, 
a single conductance step at finite bias (the transition 
from 5 = to 5 = 1), the inelastic current being three 
times larger than the elastic one. An accurate quanti- 
tative comparison with experiment of the step height is 
difficult because of the non-flat behaviour of the global 
conductance in this system (see Ref. [H); nevertheless, 
our prediction of a 300 % increase of the conductance 
at the inelastic threshold is in good agreement with the 
experimental data. For finite B, the 5=1 excited states 
split into a Zeeman structure and the present formalism 
predicts three steps in the conductance, each of them 
with a height equal to the elastic conductance. Again, 
this Zeeman limit corresponds to what is observed ex- 
perimentally for large B (see e.g. Fig. 3b in Ref. @). As 
the main feature of this system, one can stress the ex- 
tremely large contribution of the inelastic current com- 
pared to that of the elastic current in the total current: 
the present theory predicts a total inelastic current three 
times larger than the elastic one, in good agreement with 
what is observed experimentally. One can also emphasize 
that in the Co-Pc case, we only discuss the Zeeman struc- 
ture limit, so that the only ingredients in our approach 
are the value of the local spin (5 = 1/2) and of the spin 
of the electron+adsorbate system (St — 0) and all the 
predictions of relative magnitudes of elastic and inelas- 
tic currents are directly obtained from Clebsch-Gordan 
coefficients. 



V. SUMMARY AND CONCLUDING REMARKS 

The orientation of the magnetic moment of individ- 
ual adsorbatcs on a surface leads to a magnetic struc- 
ture with excitation energies in the few meV range. Re- 
cent experimental developments in the low- Temperature 
Inelastic Electron Tunnelling Spectroscopy allowed the 
direct observation of transitions among these magnetic 
states in the case of magnetic atoms adsorbed on a metal 
with a decoupling coating in between. The present pa- 
per reports a theoretical study of magnetic excitations 
induced by tunnelling electrons in metal-phthalocyanine 
molecules adsorbed on surfaces (Fe-Pc on Cu(110)(2xl)- 
O and Co-Pc stacked on Co-Pc on Pb) that have been 
the subject of recent experimental studies^. It is based 
on a theoretical framework recently introduced^ to treat 
tunnelling electron-induced magnetic transitions and on 
a DFT calculation to characterize the electronic structure 
of the Fe-Pc molecule on Cu(110)(2x l)-0 system. 



Our approach determines the strength of the mag- 
netic transitions induced by tunnelling electrons when 
the STM tip is placed on top of the magnetic atom. 
The input ingredients in our calculations are the mag- 
netic Hamiltonian describing the interaction of the ad- 
sorbatc magnetic moment with its environment and the 
total spin of the tunnelling electron-adsorbate system. 
The magnetic Hamiltonian has been taken from its pa- 
ramctcrisation using the experimental results on the en- 
ergy spectrum of magnetic levels^. In the Fe-Pc on 
Cu(110)(2xl)-O DFT study determined the elec- 

tronic structure of the adsorbed Fe-Pc molecule together 
with the symmetry and spin structure of the tunnelling 
electron. An excellent account of the experimental find- 
ings was obtained; in particular, the extremely large mag- 
netic excitation probabilities (inelastic contribution dom- 
inating over the elastic one in the tunnelling current) 
were confirmed. 

The Fe atom in the adsorbed Fe-Pc has electronic and 
magnetic structures quite different from those of a free 
Fe atom and of Fe adsorbed on CuN studied earlier^!!. 
The interaction of the d-manifold with the Pc ring and 
with the CuO substrate results in a full splitting of the 
manifold and in a spin state different from the atomic 
case. Then a Fe atom inside an adsorbed Fe-Pc molecule^ 
appears very differently in a magnetic IETS experiment 
compared to the case of an adsorbed Fe atorn^. 

The magnetic transitions appear to be much more 
probable than other inelastic processes studied earlier, 
such as vibrational excitation of the adsorbatei^— . In 
the most spectacular case (Co-Pc, see above and Ref. @), 
the inelastic contribution to the tunnelling current is 
three times larger than the elastic one. The present 
approach explains this striking difference. Since elec- 
tron tunnelling occurs on a short time scale compared 
to magnetic anisotropy, one can treat electron tunnelling 
independently of the magnetic anisotropy in the sudden 
approximation. The magnetic transitions then appear 
as the result of a change of coupling scheme for the ad- 
sorbatc spin: coupling to the adsorbate environment in 
the initial and final states and coupling to the tunnelling 
electron spin via exchange interactions for the tunnelling 
process. The strength of the magnetic transitions is then 
determined by spin coupling coefficients (such as e.g. 
Clebsch Gordan coefficients) and the present approach 
reduces to computing how a magnetism-independent tun- 
nelling current is shared among the various magnetic 
states, i.e. how a total magnetism-independent current 
is shared between elastic and inelastic parts. The impor- 
tance of a given magnetic transition is then linked to the 
weight of the initial and final states in the magnetism- 
independent collision intermediate and it can thus be 
very large. In particular, it does not depend on the 
strength of an interaction coupling initial and final states 
during electron tunnelling. Our approach is, thus, per- 
fectly well-adapted to treat situations like the present 
ones, where tunnelling appears to be dominated by in- 
clastic effects. 
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The strength of the magnetic transitions thus appears 
to be the direct consequence of the spin coupling scheme 
of the system. The variation of the magnetic transitions 
with an applied magnetic field, B, follows the variation 
of the adsorbate magnetic structure with B, basically 
the switch from a magnetic anisotropy induced by the 
adsorbate environment to a Zeeman structure, i.e. the 
decoupling of the adsorbate spin from its environment 
by the B- field action (see e.g. Fig. [6]). In this way, the 
present study of the magnetic transitions strength as a 
function of B further strengthens the knowledge of the 
adsorbate magnetic structure as it can be derived from 
the analysis of the experimental energy spectrum of the 
magnetic states. Basically, the analyses of the strength 
of the magnetic transitions and of the energy spectrum 
are probing the same properties of the system. 
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